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We develop a theory to address the non-equilibrium dynamics of phonons in a one-dimensional 
trapped ion system. We elaborate our earlier results obtained in Phys. Rev. Lett. Ill, 170406 
(2013) to chart out the mechanism of dynamics-induced cooling and entanglement generation be¬ 
tween phonons in these systems when subjected to a linear ramp protocol inducing site-specific 
tuning of on-site interactions between the phonons. We further extend these studies to non-linear 
ramps and periodic drive protocols and identify the optimal ramp protocol for minimal cooling and 
entanglement generation time. We qualitatively address the effect of noise arising out of fluctuation 
of the intensity of the laser used to generate entanglement and provide a detailed discussion of a 
realistic experimental setup which may serve as a test bed for our theory. 

PACS numbers: 


I. INTRODUCTION 

Trapped ion systems, similar to their cold atom coun¬ 
terparts, are well known emulators of specific models of 
strongly correlated quantum systems such as the Ising 
and the Bose Hubbard models Hi- A typical example 
of such a system constitutes a linear chain of ions in a 
trap The physics of the phonons, which are the 

motional quanta of these charged ions in the trap, can be 
shown to be described, under appropriate conditions, by 
the Bose-Hubbard model (BHM) Q. This leads to the 
realization of one-dimensional (ID) Bose-Hubbard model 
in such systems. The equilibrium properties and quan¬ 
tum ground states of these phonons has already been 
studied theoretically in details [H, HI. In particular, it 
was shown that such system may be the perfect test bed 
for studying the quantum phase transition between the 
localized Mott and the delocalized superfluid states of 
the phonons iii. 

The key difference between emulation of such corre¬ 
lated model Hamiltonians using trapped ions Q and ul¬ 
tracold atoms |ili3 are two fold. First, in contrast to 
their ultracold atom counterparts, trapped ion systems 
provide easier local control over the parameters of the 
model. For example, the local interaction potentials be¬ 
tween the phonons at a given site can be experimentally 
tuned for individual sites of the chain. Second, ion trap 
systems are usually prepared in form of ID chains; in 
contrast to their cold atom counterparts, it is not easy 
to prepare coherent higher dimensional configurations of 
ions in traps. The second feature makes such system 
an ideal choice of for studying correlated physics in low- 
dimension (ID) while the first aiiows for iocai controi over 
systems parameters which is very difiicuit to achieve in 
coid atoms. 

A feature of the trapped ion systems which is of key 
interest to the present study is that they aiiow one to 
study the non-equiiibrium dynamics of the modei they 


emuiate. This feature is somewhat simiiar to uitracoid 
atom emuiator; however, the trapped ion system aiiow us 
to dynamicaiiy change iocai system parameters with ar¬ 
bitrary precision. Consequentiy, one can study the non- 
equiiibrium dynamics of the emuiated modei under iocai 
protocois. Since these ion trap systems can be used to 
emuiate correiated modeis which harbors quantum phase 
transitions, they provide us with an unique opportunity 
to study the effect of iocai non-equiiibrium dynamics near 
a quantum criticai point [HI El. Moreover, trapped ions 
are ideai system to test atomic physic at high precision 
[TsI in a controiied and ciean way which makes them ideai 
units for reaiization of a quantum computer. 

In this work, we study the effect of non-equiiibrium dy¬ 
namics of the Bose-Hubbard modei, which describes the 
phonons in a iinear chain of trapped ions, under iocai pro¬ 
tocois. The kinetic energy of these phonons is described 
by a nearest neighbor hopping term whiie the interaction 
between them is iocai (on-site). Such a Hamiitonian can 
be written as 

H = J ^{blbj + h.c.) + Uihiifii - 1) (1) 

iU) * 

where J > 0 is the hopping ampiitude, Ui is the on-site 
interaction between the phonons at the site, bj is the 
annihiiation operator for the phonon at the site, and 
= blbi is the boson number operator. In what foi- 
iows we discuss the non-equiiibrium dynamics when Ui 
is ramped to —Ui on either one or two of the sites in a 
chain. We show that the first protocoi ieads to cooiing of 
the ions to their transverse motionai ground state whiie 
the second ieads to generation of an entangied Beii state 
between the two sites. We note that our proposai pro¬ 
vides an viabie aiternative for cooiing of a iong chain of 
ions which, in contrast to weii-known sympathetic cooi¬ 
ing M , does not require the presence of muitipie species 
of ions and provides a cooiing time which is indepen¬ 
dent of the eiectronic structure of the invoived ions. Fur- 
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ther, our study reveals that non-equilibrium dynamics 
can be utilized to generate pure many-body entangled 
states which are computationally relevant. 

It is well-known that the generation of entangled quan¬ 
tum states is an essential prerequisite for quantum com¬ 
putation; the present work analyzes a few ’’experimen¬ 
tally relevant” classes of possible dynamic protocols for 
generating such a state. In particular, our work here 
extends the analysis of Ref. [ij to both non-linear and 
periodic drive protocols. We show that ramp protocols 
are more efficient compared to periodic ones for cooling 
and entanglement generation and that the optimal ramp 
protocol for shortest cooling and entanglement genera¬ 
tion time corresponds to a non-linear ramp with expo¬ 
nent a ^ 0.8: Ui{t) = U^^\l — 2(t/r)^). Here is 
the initial value of the interaction at site i and r is the 
ramp time. We also show that the cooling (entangle¬ 
ment generation) can be reliably achieved by single (two) 
site addressing, z.e., the cooling (entanglement) times do 
not drastically differ if the on-site interaction on other 
sites of the chain have a near zero value. Further, we 
qualitatively address the issue of the effect of the pres¬ 
ence of noise, arising out of fluctuation of the intensity 
of the laser used to generate the on-site interaction, on 
the entanglement generation and provide an estimate of 
the maximal effective noise temperature beyond which 
such entanglement is lost. We note that single (two) site 
addressing protocols are much easier to achieve exper¬ 
imentally while the noise arising out of laser intensity 
fluctuations is an inherent part of any ion trap system. 
Thus a discussion of these issues is expected to make the 
theory developed here useful to experimentalists working 
in ion trap systems. 

The plan of the rest of the paper is as follows. First, 
in Sec. m we provide a detailed discussion of the experi¬ 
mental ion trap system which can act as test bed of our 
theoretical idea and discuss how such a system can lead 
to emulation of the Bose-Hubbard model with Hamilto¬ 
nian given by Eq. [H This is followed by Sec. Illli where 
we discuss the theoretical method used to study the non¬ 
equilibrium dynamics of the model for several local pro¬ 
tocols. In Sec. IIYI we discuss the results obtained from 
this study and identify the optimal protocols for cooling 
and entanglement generation. We provide a qualitative 
discussion of the effect noise on entanglement generation, 
summarize our results, and conclude in Sec. El 


II. ION TRAP SYSTEM 

In this section, we discuss the realization of an effective 
Bose-Hubbard model for the transverse motional modes 
in a system of trapped ion linear chain starting from their 
vibrational dynamics. For definiteness, we consider a 
string of N positively charged ions with mass m confined 
using a linear Paul trap as shown in Fig. m Sin Els. 
The RF voltage provides a confinement along the ra¬ 
dial direction whereas the axial confinement of ions is 



FIG. 1: (Color online)Schematic of the experimental setup. In 
a linear chain of 6 ions, the third ion (site no. 2) is addressed 
by an off-resonant laser to generate an on-site interaction I/ 2 . 
The axis joining the two endcap electrode is considered as the 
z-axis. 

achieved by applying dc voltage at the end electrodes. 
The Hamiltonian describing the vibrational motions in 
such a trap is given by 

N p2 

^<^ = J2^ + ^t + Vc. ( 2 ) 

i=l 

Here is the momenta of the ion, Vc is the coulomb 
potential between any pair of the ions which can be ex¬ 
pressed as 

1 ^ 2 

vc = - y ^ 

2 \J(Xi - Xj)‘^ + ivi - VjY + {Zi - ZjY ’ 

( 3 ) 

where i/i, and Zi denote the coordinate of the ion 
in three orthogonal directions, and Vt is the trapping 
potential, which determines the equilibrium positions of 
the ions and is given by 

i=l 

Here are the radial and longitudinal trap fre¬ 

quencies. In what follows, we are going to work in the 
regime oJx^^y ^ so that the ions are strongly confined 
in the radial direction. In our setup, the chain direction 
is taken to be along 2 : as shown in the Fig. [TJ The oscil¬ 
lations of ion along the radial directions, in this regime, 
is small compared to the inter-ion distance. 

Both the trapping potential and Coulomb repulsion de¬ 
termines the equilibrium position of ions. In the regime 
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described above, the Coulomb potential is much weaker 
compared to the radial confinement, and the displace¬ 
ments of ions in any direction can be considered to be 
decoupled from the other two orthogonal directions al¬ 
lowing us to discuss the vibrational modes in the radial 
direction only. To this end, we expand Vc as a function of 
radial displacement around the equilibrium ion position 
and obtain the Hamiltonian 


p2 

U \^i,x , ^ 2 e 2 

i=l 


e^{5xi 




14 p 

(*>i) 


( 5 ) 


where \d^- \ = \z^ — z^\ = dQ\iz — jz\ is the mean distance 
between and ion along the 2 ;—direction, iz{jz) are 
the site number of the site, do is the inter-site 

distance between neighboring ions, and Sxi is the radial 
displacements of the ion around its radial equilibrium 
positions. The second quantized form of this Hamilto¬ 
nian can be obtained following standard methods 0 
and takes the following form {h = 1): 


N N 

Hx = jjj + h.C.) + + ^x,i)Y^i^i (6) 

i,7 = l i=l 

i>j 

Here bi denotes the annihilation operator for phonons 
(quanta of vibrational displacement modes) in the ra¬ 
dial direction, and the hopping amplitude Jij is induced 
by the Coulomb interaction. The effective trapping fre¬ 
quency, ujx,i depends on the ion’s positions, as well as the 
tunneling amplitudes Jij and is given by 


(j^i) 


2muJx\cPiX'^'^^ 


In this context, it is customary to define the dimension¬ 
less parameter px as the ratio of the Coulomb energy 
between near-neighbor sites to the trapping potential en¬ 
ergy: 


Px = e^/(mujldl). ( 8 ) 

The parameter px characterizes the normal modes of the 
system in radial direction; when px » 1, the Coulomb 
repulsion dominates over the trapping potential. In this 
case, the long-range nature of the Coulomb interaction 
becomes important and the effect of the collective mo¬ 
tion along the axial direction is expected to play impor¬ 
tant role in the collective motion of the ions and can not 
be ignored. However, in the limit Px « 1, the Coulomb 
interaction is effective between nearest ions only. In this 
regime, which we focus on, the number non-conserving 
term in the effective Bose-Hubbard Hamiltonian such as 
(pihj^ b\h'^j) can be neglected since they lead to an en¬ 
ergy cost which is larger than the Coulomb energy. We 


have used this fact to arrive at Eq. [6l Also, since the 
Coulomb interaction is effective only between the near¬ 
est neighboring sites, one can express Ji^j in terms of Px 
as 

Jij nZ Jdi^j±^\^ J = Px^xl‘^' (9) 

Note that for /3a, ^ 1, the effect of on-site phonon- 
phonon interaction is negligible due to the small ampli¬ 
tude of anharmonicity in the trap which arises from the 
imperfection of trap architecture. However, in a realistic 
system, this effect can be enhanced and controlled by in¬ 
troducing off resonant standing wave along the transverse 
direction of the trap axis. The bosonic Hamiltonian, in 
the presence of such a standing wave, has an additional 
term 

Hi = fY^ cos^{ kxi + ^S), (10) 

where F is the peak AC stark shift and k is the wavevec- 
tor of the wave. This controllable anharmonicity brings 
the phonons closer by making unequal the otherwise 
equally spaced energy levels of radial trap potential. To 
quantify this effect, one may define the dimensionless 
Lamb-Dicke parameter r]x = ky^h/2mujx « 1; in terms 
of ? 7 , it is possible to write Hj as 

Hi = FY,[^Wx{bl + bif + ^vM + bi)^+HVxn (H) 

i=l 

In the regime where Frj^ << cjx, the effective interaction 
between the bosons is quartic and can be written as 

N 

Hi = {-l)^F4/3j2ibi + blf ^ ^ Uib]^bl (12) 

i=l i 

where in arriving at the last relation we have used the 

9 F 

fact that the number non-conserving terms such as 6^,6] 
can be neglected since they corresponds to much higher 
energy scales as discussed earlier. Also, it can be shown 
that the rest of the phonon conserving terms will induce a 
small shift in the radial potential which amounts to a re¬ 
definition of the the overall potential. The above discus¬ 
sion leads to the final form of the Bose-Hubbard Hamil¬ 
tonian given by Eq. [T]with J > 0 and Ui = 2{—lYFr]Y 
The on-site interaction strongly depends on the strength 
of the standing wave and relative posit ion (J = 0 or 1) 
of ion on the standing wave. Note that it is possible to 
have Ui to be repulsive or attractive which is one of the 
key facts used in the present work. Einally, it is useful 
to note that the sign of the hopping term in the present 
Bose-Hubbard model is positive in contrast to its ultra¬ 
cold atom counterpart where such a term comes with a 
negative sign. 


III. METHOD FOR STUDYING 
NON-EQUILIBRIUM DYNAMICS 

In this section, we outline our method for study¬ 
ing non-equilibrium dynamics of the ID Bose-Hubbard 
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model (Eq. [T]) derived earlier. We begin by specifying 
the protocols for the dynamics that we study. We start 
with the system in the ground state of the Hamiltonian 
given by Eq. [T]for a fixed J and Ui. We consider a chain 
of N ions; the choice of N depends on experimental set¬ 
ting which is to be discussed in Sec. [Vl 

The dynamics is induced in the model by changing the 
local interaction parameter to on a 

one or more selected site(s). Eor cooling, we change Ui at 
any one of the N sites, while for generating entanglement 
this change is made for local interactions on two of the 
sites of the chain. The protocol used to make such a 
change are either linear or non-linear ramp or periodic 
and are given by 


with the boundary condition Cck(O) = To solve these 
equations, it is convenient to rewrite 

H{t) = i7(r) +Ai7(t), 

Aff(t) = Y,{Ui{t)-Ui{T))hi{fn-l). (17) 

i 

With this choice, one obtains the final set of equations 
for Ca{t) to be 

{ihdt - Ea)Ca{t) = 

= {P\AH{t)\a). (18) 


Ui{t) = — 2{t/r)^], 0 < t < r, (ramp) 

= cos(7rt/r), 0 < t < r, (periodic)(13) 

where the exponent a specifies the power-law which is 
followed by the ramp and a = 1 indicates linear ramp 
protocol studied in Ref. [TH. We note here that r ^ 0 
indicates the sudden quench limit; such dynamics for the 
model was studied in Ref. In case of periodic ramp, 
we would be using T and cjq to denote the period and 
frequency of the ramp which is related torasr = T/2 = 
tt/cjo- 

To study the dynamics, we begin at t = 0 assum¬ 
ing that the system is in the ground state of the Bose- 
Hubbard model given by Eq. [T]with > 0 at 

all sites and a fixed at each site. Since we work 

in a regime where changing the number of phonons is a 
high-energy process, we shall study the dynamics at fixed 
number Nq total phonons (bosons). Thus the Hilbert 
space of the bosons can be truncated to keep those states 
for which Ui < Nq for any site i. We use exact diago- 
nalization method for the finite-size system within this 
truncated Hilbert space to obtain the energy eigenstates 
\a) and eigenvalues for H{t = r). This amounts to 
a choice of basis; in terms of these eigenstates, one can 
express the initial ground state (at t = 0) 

= (14) 

a 

where the coefficients denote the overlap of the initial 
ground state of the system, also obtained using exact 
diagonalization in the same truncated Hilbert space with 
H = H{t = 0), with \a). 

The time-dependent Schrodinger equation for the sys¬ 
tem wavefunction 

= (15) 

a 

governing the dynamics of the system now reduces to 
equations for time evolution of Ca{t): 

ihdt Y,Ca{t)\a) = H{t)J2ca{t)\a) (16) 


The set of coupled equations for Ca{t) is then solved nu¬ 
merically leading to an exact numerical solution for the 
time-dependent boson wavefunction |'0(t)). We note here 
that it is clear from the structure of Eq.[T8]that Ca{t) can 
be written as 

Ca{t) = 

ihdtCait) = (19) 

13 

Having obtained the wavefunction |'0(t)), one can com¬ 
pute the expectation value and equal-time correlation 
functions corresponding to any boson operator as a func¬ 
tion of time. This is given, for a local boson operator Oj, 
as 


Oj{t) = 

al3 

c;r^O(i) = (20) 

a/3 


where n, m are arbitrary integers and j, k denote site in¬ 
dices. Although our method is powerful enough to al¬ 
low us to compute the expectation and correlation corre¬ 
sponding to any bosonic operator O, in what follows, we 
shall be mainly interested in obtaining these quantities 
pertaining the boson annihilation operators bj and the 
number operators hj = b^jbj. 

Before ending this section, we briefly comment about 
the single/two-site addressing protocol which we shall use 
in Sec. [IVl Eor such protocol at t = 0 only 

on the site(s) where they are dynamically changed dur¬ 
ing the protocol. Eor the rest of the sites Ui = 0 at all 
times. This protocol has the advantage of experimen¬ 
tal simplicity, and as we shall, see leads to faster cooling 
and entanglement generation. We note here that the for¬ 
malism developed here can be applied for this protocol 
without any major modification. 
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IV. RESULTS 

In this section, we present a detailed account of the re¬ 
sults obtained by numerical solution of Eqs. [19] followed 
by evaluation of appropriate expectation value or cor¬ 
relation function by computing Eq. |20l In Sec. IIV Ai 
we discuss these results in the context of cooling. This 
is followed by Sec. IIVBI where we discuss entanglement 
generation. 


A. Cooling of ions 

To achieve cooliM of the ions, we follow the protocol 
developed in Ref. [12 and change the on-site interaction 
on one of the lattice sites in the chain from 

to with a fixed ramp rate r In Ref. 

E we have studied the linear ramp protocol. Here we 
concentrate on non-linear ramp protocol characterized by 
an exponent a (Eq. [T3|) with the aim to optimizing the 
cooling time. Eor the sake of definiteness, we shall choose 
a chain of L = 8 site and choose to ramp the interaction 
on the second site of the chain {i = 2 ). We note that this 
choice is made for the sake of definiteness and do not lead 
to loss of generality. In what follows, we study the time 
evolution of the phonon population on the second site: 
N 2 /N = ( 62 ^ 2 )/^- In the rest of this section, we shall 
set the total phonon number N = 4. 

The results of such a study are shown in Eig. |2l. Eig. 
O In Eig. [21 we plot the phonon number N 2 at the end 
of the ramp (t = r) as a function of the ramp time r for 



FIG. 2: (Color online) Expectation values of N 2 /N as a func¬ 
tion of the ramp time r and exponent a. Linear ramp is repre¬ 
sented by a = 1 while the other three are for exponents 0.1, 0.8 
and 3 respectively. For all plots the initial =0.2. The 

maximal phonon transfer occurs when a = 0.8(red). 



FIG. 3: (Color online) Comparative plot of cooling in case 
of all site except site no. 2 applied with a repulsive on-site 
interaction {Ui > 0) and no on-site interaction {Ui = 0). The 
dashed-dot curve (black) is same as in Fig. [2] for the optimal 
exponent a = 0.8 while the solid-line curve is for no on-site 
interaction except at site number 2. The cooling time remains 
similar while the final state is cooler in case of no interation 
as comapred to interaction applied to all sites. 


several a and = 0 . 2 . We find that, in accordance 

with our expectation, N 2 approaches N for slower ramps 
indicating the migration of the phonons to the second 
site of the chain leading to cooling of the other sites. 
By studying the evolution of the phonon number on the 
second site, we find that N 2 approaches N in the shortest 
possible time for a 0.8. Thus our study reveals that 
a non-linear ramp with a = 0.8 is the optimal protocol 
for cooling. We find that with = 2 . 8 kHz on all 

sites, the linear ramp achieves N 2 /N = 0.9(0.97) for t = 
12(26)ms; in contrast a ramp with a = 0.8 with same 
leads to N 2 /N = 0.9(0.97) for t = ll(24)ms. 

Next, to study other aspects of the cooling dynam¬ 
ics, we resort to single site addressing, z.e, we keep the 
initial interaction parameter = U^^^ 6 i 2 finite only 

at the site where it is to be dynamically changed to 
interaction is set to zero for phonons on rest of 
the sites. The motivation behind such single site ad¬ 
dressing is that it is experimentally a lot simpler; further 
as shown in Fig. [3l such single site addressing does not 
lead to appreciable change in the nature of N 2 /N as a 
function of r. In fact, the cooling time corresponding 
to N 2 /N = 0.97 can be further reduced by such single 
site addressing: for a = 0.8 and = 2.8 kHz, one 

reaches N 2 /N = 0.9(0.97) in t = 10(19)ms. In the rest 
of this section, we shall thus concentrate on single site 
addressing. 

We further look for optimization of the parameter 
within single site addressing. To this end, we 
study the behavior of N 2 /N at the end of the ramp 
{t = r) as a function of for several representa- 
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tive values of r and for a = 0.8. Our study reveals that 
the maximal value of A^ 2 /^ occurs at ^ 0.2. This 

leads us to identify r > 500, a 0.8, and 0.2 

as optimal ramp parameters for achieving cooling. 



FIG. 4: (Color online) Fractional photon trans¬ 

fer N 2 /N as a function of initial BH parame¬ 
ter for three different ramp time r = 

100 (square,black), 300 (triangle, blue), 500 (circle, red) 
(in units of 1/\U^^^\). The ramps are non-linear with 
optimized value of o = 0.8. In all these cases, the on-site 
interaction has been applied to a single site 2. 


with a near-zero slope, followed by a steeper increase at 
intermediate and later times which leads to worse over¬ 
lap compared to a ramp with optimal exponent.We also 
note that the ramp protocol advocated here is not nec¬ 
essarily the absolute optimal protocol for the problem 
0 ; however, it is certainly the most efficient among the 
experimentally accessible protocols discussed here. 

B. Entanglement generation 

In this subsection, we discuss dynamical entanglement 
generation using the ramp protocol. In doing so, we 
choose to change the sign of the interaction following 
the ramp protocol given in Eq. [13] on two symmetrically 
placed sites on the chain. For definiteness, and with¬ 
out loss of generality, in the rest of this section, we shall 
choose a chain of length L = 6 with site index 0 < i < 5 
and change the on-site interaction for k = 1 and 1 = 4. 
We measure the cross correlation between the these two 
sites which is given by 

Cki{t) = {i^{t)\{blbi)^\'tp{t))/N\ (21) 

We note that if we carried out the protocol of changing 
the sign of interaction for sites k and I adiabatically (with 
infinitely slow rate; r ^ 00 ), Ci 4 (t/ = 00 ) = = 

1/2. Such a correlation, in the present context, pertains 
to the Bell state between the sites k and I given by 


Finally, we compare ramp protocols to periodic proto¬ 
cols, described by Eq.[T3l for identifying the optimal pro¬ 
tocol for cooling. In what follows, we study the evolution 
of A^ 2 /^ under such a protocol for =0.2 and with 

single site addressing. The result of such a study is shown 
in Eig. |5l In the top left panel of Eig. |5l we plot N 2 /N at 
the end of the drive as a function of the drive frequency 
ujQ. We find that N 2 /N has an oscillatory behavior as 
a function of uoq for small changes in ujq; however, for 
a larger variation of cjq, h decreases with increasing ujq. 
The rest of the panels of Eig. [5] shows the time evolution 
of N 2 /N for a wide range of c^q. Erom these plots, we 
find that one needs a much longer cooling time if peri¬ 
odic protocols are used; for example, N 2 /N = 0.9(0.97) 
for t = 81(88)ms with uoo = 0.0012 in terms of 

= 2.8kHz and singe site addressing. Thus we con¬ 
clude that ramp protocols perform better than their peri¬ 
odic counterparts for the purpose of cooling. The origin is 
this behavior can be qualitatively understood as follows. 
Eor efficient cooling, one needs to maximize the overlap 
of the final system wavefunction with the ground state 
wavefunction. Such an overlap depends on the slope of 
Ui{t) at each instant of time. This can be understood 
by noting that a near zero slope of Ui{t) is expected to 
reproduce the adiabatic evolution which will lead to max¬ 
imal overlap. The slope of Ui{t) can be optimized more 
efficiently for ramp protocols by varying the exponent 
a. Eor periodic protocol such as the one described by 
Eq. [131 fhe short time {ujot 1) evolution always occur 


l^Bew) = ^(|0Ar0000) + |0000Ar0)), (22) 

V2 

since it yields C 14 = Thus in the adiabatic limit, 

1 7 /) generated by the protocol described above is expected 
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FIG. 5: (Golor online) Plot of N 2 /N for periodic ramps: Ini¬ 
tial = 0.2 and single site addressing done for all plots. 

Top-left panel is N 2 /N vs loq and rest of the figures show 
N 2 /N as a function of t for three different values of Ct;o : 0.1, 
0.01, and 0.001 respectively. 
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to have perfect overlap with the Bell state. However, if 
the ramp is done fast, as shown in Ref. in Ci 4 {tf) does 
not reach even when tf ^ r. Since one of the 

motivations behind generation of entangled state is per¬ 
forming quantum gate operations in the shortest possible 
time, it is therefore useful to find optimal ramp time r 
and power a for which Ci 4 {tf) reaches a significant frac¬ 
tion of This issue has been investigated for linear 

ramp protocols in Ref. [Ill; here we generalize such study 
to non-linear ramp protocols. 

In attempting such generalization, we find that in ana¬ 
log to the dynamic protocol for cooling of ions, the en¬ 
tanglement generation protocol also yields the shortest 
time for (a = 0.8. So in what follows, we concentrate on 
a non-linear ramp with a = 0.8. Furthermore, in analogy 
to the cooling protocol, we resort to two-site addressing, 
z.e., we keep a finite initial interaction only the on two 
sites where we change the interaction parameter; for the 
rest of the sites, = 0 at all times. We note that, as 
shown in Fig.[6l the two-site addressing do not change the 
dynamics of Ci 4 {t) in any appreciable manner; in fact, 
such a protocol leads to faster plateauing of Ci 4 {t) which 
signifies a faster entanglement generation. Thus we shall 
restrict ourselves to the two-site addressing protocol for 
the rest of this section. 

Next, we study the maximum attained value of the cor¬ 
relation, in a given ramp with exponent a = 0.8 

as a function of the ramp time r and J/U^^\ The corre¬ 
sponding plot, shown in Fig.[3 indicates that it is indeed 
possible to reach very close to the maximal value 1 /2 of 



FIG. 6: (Color online) Phonon-phonon correlation between 
sites 1 and 4 Ci 4 as a function of the ramp time r in units of 
for a non-linear ramp of exponent a = 0.8. The dashed 
(red) curve is for = 0.15 applied to sites 1 and 4 only, 

while at the other sites Ui — 0. The solid-line (black) curve 
is for all sites initially at =0.15 while ramp is applied 

to only sites 1 and 4 while others are kept 3.1 J/U = 0.15. As 
in case of cooling (see Fig. 0, application of on-site interac¬ 
tion to all sites during the ramp doesnot make any significant 
difference. 



FIG. 7: (Color online) Entanglement generation in terms of 
cross-correlation as a function of ramp time and initial 
value for non-linear ramp of exponent a = 0.8 and two site 
addressing only. 



FIG. 8: (Color online) The real part of the density matrix for 
six site linear chain. In this case the correlation Cki{0) is for 
= 0.15 applied to only sites 1 and 4 which has been 
considered as the initial state. 


Ci 4 within a finite time; the optimal parameters for this 
turns out to be r = b00h/\U^^^\ and = 0.15. To 

show that such a correlation indeed represent realization 
of a Bell state, we compute the cross correlation between 
all pairs of sites, Cki for r = b00h/\U^^^\, a = 0.8 and 
= 0.15 for sites i = 1 and i = 4. At t = 0, before 
the ramp is carried out, we find from Fig. [8l that Cki{^) 
is finite and has comparable magnitude for most pair of 
sites. This situation is to be contrasted with the behav- 
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ior of Cki{t = r); as shown in Fig. [9l Cki{t = r) is close 
of to 1/2 for k = 1,4 and / = 1,4 and is close to zero 
otherwise. This is clearly a signature of generation of 
Bell state through dynamics. We note that this method 
may also be used to generate entangled states involving 
multiple sites such as the Greenberger-Horne-Zeilinger 
state HI; however generation and maintaining stability 
of these states are expected to be much more complicated 
experimentally. 



FIG. 9: (Color online) Tomograph for J/\U^^^\ = 0.15 with 
< 0 for /c = 1 and 4 and > 0 otherwise. The other 
parameters are T\U^^^\/h = 500 and a = 0.8. 

Before ending this section, we would like to make a 
comment regarding the requirement that the sites be¬ 
tween which such cross correlation is generated need to 
be symmetrically placed for a finite chain. This require¬ 
ment stems form the fact that the on-site energies of 
these two sites need to be identical for generating a Bell¬ 
like state. These on-site energies have contribution from 
virtual boson hopping processes which are or 

higher. Since the boson hopping processes are cut off 
by chain ends, these contributions vary with the posi¬ 
tion of a site; thus two sites which are asymmetrically 
placed with respect to the chain ends would have differ¬ 
ent on-site energies due to different contributions from 
virtual boson hopping processes. Since for creating en¬ 
tangled states involving bosons in two different sites, one 
needs to have equal on-site energies of the two sites to 
a high-degree of accuracy, these sites need to be placed 
symmetrically with respect to the chain ends. Indeed, as 
shown in Fig.[T0l a variation of on-site energy induced by 
changing the ratio from one, leads to a decay 

of the entanglement. This inherent fragile structure of 
the entangled state motivates us to at least qualitatively 
discuss the effect of the presence of noise (which is an 
integral part of any experiment carried in ion trap sys¬ 


tems) on entanglement generation. We shall discuss this 
in the next section. 


V. DISCUSSION 

In this work, we consider cooling and entanglement 
generation for an ion trap system with a chain of L sites 
with N phonons. The numbers L and N can be con¬ 
trolled and they vary within the range 2 < L < 14 and 
0 < A/" < 10 in typical experimental setups. For cooling, 
we choose any one (i^^) the L sites with at 

t = 0 and vary Ui{t) according to a ramp protocol so 
that at the end of the ramp. This leads 

to migration of phonons to the sites leading to cool¬ 
ing of the rest of the sites in the system. We find that 
Ni/N ^ 0.97 for r ^ 19ms for a typical \ = 2.8kHz 
and ^ 0.15. Such a cooling time pertains to a 

chain of L = 8; one expects the cooling time to increase 
linearly with the chain length since the migration time of 
the phonons ^ LhjJ increases linearly with this length 
for a fixed J. 

Next, we come to the issue of entanglement genera¬ 
tion. As we noted in Sec. lIVBi the dynamical genera¬ 
tion of entanglement crucially depends on the equality 
of the on-site energies of the two chosen sites i and j of 
the chain; Cij falls down rapidly from its maximal value 
of 0.97Cf7^^ when these on-site interactions differ from 
each other. Thus the entangled state generated during 
dynamics is expected to be fragile against any relative 
fluctuation of the on-site interactions Ui and Uj. Such 
fluctuations are integral part of a realistic experimental 
setup since they originate from the fluctuation of the in¬ 
tensity and/or frequency of the lasers used to generate 
the on-site interactions. Thus any realistic protocol for 
dynamic generation of entangled state needs to address 
the effect of noise on Cij; in what follows, we qualitatively 



FIG. 10: (Color online) Plot of Ci 4 (t/) with the ratio of two 
on-site interaction strengths Ui and Ua reflecting the sensitive 
dependence of the entanglement on the on-site energy. 
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address this issue. 

Assuming that slow drifts can be controlled by suitable 
locking, the time scale of the laser fluctuations ^ 100ns 
are much shorter than typical entanglement generation 
time r ^ 20ms. Thus in any typical measurement, one 
expects the effect of noise to self-average; an estimate 
of Cij can thus be obtained from its average value over 
several noise realization. A direct estimate of Cij thus re¬ 
quires solution of stochastic differential equations given 
by Eg. [T8l with Ui{t) Ui{t) + 6 Ui{t), where 6 Ui{t) de¬ 
notes the random fluctuation of the on-site interaction 
on the site. We do not attempt to solve this problem 
here. Instead, to obtain a qualitative understanding of 
the effect of noise, we consider a two site problem with 
two phonons (L = 2, = 2) which are described by the 

Hamiltonian 

Hi{t) = J{b\b 2 + h.c) + U{t)[niihi-l) 

+{l + SU/U{t))n 2 {n 2 -l)]. (23) 

Here 61(62) are the phonon annihilation operators on 
sites 1(2), hi = 6|6i is the phonon number operator, 
J is the hopping amplitude between the sites, U{t) = 
[/(o)(i _ 2(^/i-)«) is the on-site interaction which is var¬ 
ied from to within a ramp time r, and SU 

is a random number representing the relative fluctuation 
between the on-site interactions at site 1 and 2. In what 
follows we shall choose 6 U to have a Gaussian distribu¬ 
tion exp[—where the variance r]o = denotes 

the effective noise temperature of the system. 

The Hilbert space of the above model can be charted 
out in terms of the number states of bosons on each site. 
For simplicity we are assuming that the ions start in the 
entangled state and the noise is imposed at the two sites 
of entanglement. For N = 2^ the Hilbert space consists 
of three states which are given by |0) = |ni = l,n2 = 1), 
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FIG. 11: (Color online) Plot of Ci 4 (t/ = |) with 

the effective noise temperature 770 in units of Ho showing the 
gradual decay of entanglement with larger noise temperature. 
The initial value of J/Uq is chosen to be 0.2 for this plot. 


|1) = \ni = 2,77-2 = 0), and |2) = |77-i = 0,77-2 = 2). 
Thus the state for the bosons at any time t during the 
dynamics and for a given realization of the noise SU can 
be written as 

2 

= (24) 

a=0 

where the equations for the coefficients Ca{t) can be 
obtained from the Schrodinger equation ihdt\'ip{t)) = 
and are given by 

ihdtco = \/2J(ci+C2) 

ihdtci = U{t)ciV 2 JC 0 

ihdtC 2 = [U{t) ^ SU]c 2 ^ V 2 JC 0 . (25) 

Eqs. [25] can be solved numerically to obtain Ca{t) for a 
given disorder realization. The initial conditions for solv¬ 
ing such equations can be obtained from finding their 
value for the system ground state at a fixed value of 
J/U^^\ This can be done by numerical minimization 
of the system energy E = ('0(O)|iJ(O)|'0(O)). 

Having obtained (t) for a given disorder realization, 
one can compute the cross correlation C 12 given by 

Cuit) = {'>Pit)\{blbif\i’{t))/2 = dlc2 (26) 

One then averages C 12 over several disorder realization 
and obtain Ci 2 {t) = 2Re[(Cf2(^))d], where the average 
is taken over M = 100 realization of disorder. The be¬ 
havior of C12 (t) as a function of the effective noise tem¬ 
perature 770 is shown in Fig. [TT] for Ji/U^^^ = 0.2 and 
t = 5006-/1I. We note that the entanglement gradu¬ 
ally decays from its value 0.9 for 770 = 0 to around 0.2 for 
770 = with increasing 770. In a typical experimental 
setup, the variation of the laser intensity and/or phase 
is around 1% and this leads to around 3% fluctuation in 
the value of U. This suggest that the experimental noise 
might lead to a reduction of ^ 10% value of the cross 
correlation as obtained from the simple naive model. A 
more detailed analysis of this phenomenon which is ex¬ 
pected to provide an accurate estimate of the entangle¬ 
ment reduction in a more realistic situation is left for 
future work. 

In conclusion, we have studied the non-equilibrium dy¬ 
namics of phonons in an ion trap and have shown that 
such dynamics, for carefully chosen protocol, may lead 
to both cooling and entanglement generation. We have 
provided a detailed analysis of theoretical method used 
to study such dynamics. Using this, we have identified 
the non-linear ramp protocol with exponent a = 0.8 to 
be optimal for both cooling and entanglement genera¬ 
tion. Finally we have provided an estimate of cooling 
and entanglement generation times based on our analy¬ 
sis which may be confirmed in experiments and have also 
qualitatively discussed the effect of the presence of laser 
intensity fluctuations on the entanglement generation. 
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